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ABSTRACT 
The   statistical     model  assumes   a   statistical   distribu- 
tion of   nuclear   density   levels   in  the   decay   of   the   compound 
nucleus.      The   densities   of  nuclear   levels    can  be   described 
by  a   statistical  formula   in  a   region   containing    many   energy 
levels.      This   formula   can  be   used  to   predict   the   cross- 
sections   for   photonuclear   reactions   within  the   energy 
range  0    to   30   MEV .      Provided  assumptions   are   made 
neglecting  the   coulomb    effect,    spin,    parity,    angular    mom- 
entum,   and  tertiary  and  higher   order   reactions,    these 
predictions    can  be   made  by  a   digital   computer   and   the 
results   plotted   on  a   graph   with   experimental   data   for 
comparison.      The   digital    computer   program   is   used  to   find 
the   statistical   level   density  formula   that   most   closely 
approximates   the   cross-sections   obtained   experimentally. 
The  formulas   used  thus   far   do   not   approximate   the   exper- 
imental  data   closely   enough.      Some   formulas   for   the 
density  functions   give  a   shape  for   the   cross-section  as 
a  function  of   energy   that   corresponds    well   with   the 
experimental    curve,    but   the   ordinates   are   too    small.      The 
formula   P(E)    =   C    E    EXP    (A   En)    is   the  best   one   found 
so   far . 
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1.  Introduction 

A  description  of   a  nuclear   reaction  in  terms   of   the 
decay   of   long— lived   equilibrium   systems   which  neglect   phase 
relations   can  be   called  a    "statistical    model".      It   is 
recognized   that   statistical  theory   qualitatively  predicts 
the   major   portion  of   nuclear   reaction   cross    sections   and 
in  this   thesis   a   comparison  is   made  between   the   statistical 
model  and   experiment.      The  reason  for   this   is   primarily 
that   the   decay  of  an   equilibrium   system   is   governed  by 
the  available  phase  space.  1      The  use   of   nuclear   models   is 
an  attempt   to    explain  in   tangible  form   the   observed 
phenomena   characteristic   of  nuclear   reactions. 

To    explain  large  resonances   in  neutron      cross-sections, 
Niels   Bohr   in  1936   introduced  the   concept   of   a   compound 
nucleus   which  is   a   many-body   system    of   strongly  interact- 
ing particles   formed  by   the   combination  of   an  incident 
particle  a.  with  a   target   nucleus: 

X   +   a— ►  C* 
In   such  a   system    the  incident   particle  a  has   a   short   mean 
free  path  and   shares   its    energy   with  the   other   particles 
of    the   system    C*    so    that   it   cannot   be  re-emitted  until, 

Ericson,    Torleif ,    The   Statistical  Model  and   Nuclear 
Level   Densities.      Advances   in   Physics,    v.    9>    no.    36, 
Oct.    I960:    i+30. 


as   a   result   of   further   exchanges,    sufficient    energy   can 

again  be   concentrated  on   this   or   on  a   similar   particle. 

Should  the  incident   particle  be  a   slow   neutron,    this    could 

be  a   length   of   time   sufficient  to   permit   the   emission  of 

radiation.      The   capture  process   is   complete,    and  the 

compound  nucleus   C   is   formed  in  its   ground   state   either 

by   emission  of  a  single  photon 

C*  — ►  C   +   V 

or  by  the  emission  of  a   sequence  of  quanta  in  cascade. 

This   explains  the  predominance  of   the  capture  reaction 

and  the  suppression  of   elastic   scattering. ^ 

The  progress   of   many  types  of  nuclear  reactions 

can  be  described  in  terms  of   the  compound  nucleus  as  a 

two-stage  process: 

X   +  a-^C*— ►Y   +  b 
in  contrast   to   a  single-stage   process: 

X  +  a  — ►  Y   +  b     3 

Niels    Bohr   observed   the  usefulness   of   dividing  a 
nuclear  reaction  into   two   states: 

(1)  the  formation  of   the   compound   system 

(2)  the    disintegration  of   the   compound   system 
into  the  reaction  products. 

2Blatt,    <John  M.,    and  Weiskopf,    Victor   F., 
Theoretical    Nuclear   Physics.      «John  Wiley  and   Sons,    New 
York.    1952:    528-529. 

3Blatt  and  Weiskopf:    528-529. 
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These  two   states   can  be  treated  as  independent  processes 
in  the   sense   that   the   mode  of   disintegration  of   the   com- 
pound system    depends    only   on  its    energy,    angular   momentum, 
and  parity,    but  not  on  the  specific   way  in  which  it  was 
formed.      This   assumption  has   proved   to   be  valid   or   approx- 
imately valid  for   many  cases. ^ 
2.         Theory   of    Nuclear   Reactions 

a .      Channels 

Consider   the  following  nuclear   reactions: 


a      (elastic   scattering) 
a 
Y  +  b 
yS   +   c'        etc. 


All   of   these  reactions   except  that   of   elastic   scattering 
can  be  subdivided  according  to   the  quantum   state  of   the 
residual  nucleus  and  the   emerging  particle.      One  can 

denote  the  states   of  these  nuclei  by     o(    j  {?     .   Y     , , 

and  the  states   of   the  incident   or   emerging  particles  by 

d\      ,  {9     ,    V        , ,.      Should  particles  a,b,c, ,    be 

elementary,    states     o(    .    /? "       'y"i 1    refer   to   their 

spin  orientation,    and  the  following   reactions   result: 


V +  xdki 

Y^i  +     br,» 


^Blatt  and  Weiskopf :    341. 


0(     and  c^"   are  states    of   the  target   nucleus   and  the 
incident   particle;    thus,/9l   and  ,#",    or  y' and    'y"    can   denote 
any  quantum   state  of    X   and  a   or   Y   and  b   respectively 
which  can  occur  in  this   reaction.      The  conservation  laws 
of    energy,    angular   momentum,    and  parity   restrict   the 
possible  combinations  of  A  '  and#  ",  -y"1  and  y"      ,    etc. 
Any   such  possible   combination  of   residual  nucleus   and 
emerging  particle,    each  of   which  is   in  a   definite   quantum 
state,    is   called  a   reaction  channel.    These   channels   are 

denoted  by   single  Greek  letters  c{t /9  »  Y~i »    which 

consist   of  both    indices       o^'   and  </Jr   or/P'     and^"'.    Channel 
X  ,i  +  a  .  "  is   called  the  entrance   channel  or  initiating 
channel  of   the  above  reaction.^ 
b.      Reaction  Probability 

According  to   the  Bohr  assumption,    one   can   write 
the   cross-section  of   a  nuclear   reaction  X(a,b)Y   in  the 
form: 

6la,b)    =Q(a)     Qc  (b) 
where     6^(a)    is  the   cross-section  for  the  formation  of  a 
compound   system   by  particle  a  incident   upon  the   target 
nucleus   X.    £^(b)    is   the  probability  that   the   compound 


5Blatt  and  Weiskopf:    313 . 


system   C,    once  formed,    decays  by  emission  of  a  particle 

b,    leaving  a  residual  nucleus  Y.   £-(b)    is  a  pure  number; 

the   compound   system    C    must   eventually   decay  in  some 

special   way;    G-Q{b)    is   the  probability  for   this    special   mode 

of  decay." 

It  is  useful  to   describe  the  reaction  X(a,b)Y  in  more 

detail  by   considering  the   cross-section    S'C^j^J  corresponding 
to  a   specific   entrance  channel  cJ\  and  a  specific   exit 
channel  j9  •      Therefore,    one  specifies   the  quantum   states 
of  all  reaction  pairs  before  and  after  the  reaction.      Thus: 

G(vL,p)  =   6TC  (d)  6c  0)  where  6Tc(d)     is 

the   cross-section  for  the   formation  of   C   through   channel 
q^    ,    and  &c ifi)    is   the  probability   of   C   decaying   through 
channel  B  . 

The  probability  Q-{B)    that   the   compound  nucleus    C 
will   disintegrate  through   a   specific    channel  jS  can  be   written 

in  the  form 

>3l 


<&)  ~-  IC  ^M) 


where  k  =  X  and  the  sum  is   extended   over  all   channels  V 
into    which  the   compound  nucleus    can   decay.      Thus   the 


Blatt  and  Weiskopf:    342. 
7Blatt  and  Weiskopf:    3^2. 


cross-section  for   the  reaction  dk- ►^  can  be   calculated  if 
the   cross-sections  6!  (Vj for   the  formation  of   the   compound 
system   by  all  possible   channels   are  known  and  if   the  Bohr 
assumption  is   valid.      The  probability  &(&) '1S   energy   dependent, 
and  its  value  for  a   compound  nucleus   excited  to  an  energy 
Ec    can  be   obtained  from   the  above   expression  provided 
the  cross-sections  6^(y)  and  the  wave  numbers  kr     are 
known  for   those   channel   energies  Sy     which  produce   the 

o 

compound  nucleus   with  an   energy   Ec« 

c.      Decay  of   the  Compound  Nucleus 

The  kinetic   energy  in   channel  f9  is   given  by 
C   -Z.C     +    /Q        where  Q      is   the    "Q"   value  from   channel 
d^    to   channel  j9    .      Because  most  of   this   energy  is   the 
kinetic   energy   of  the   emitted  particle  b  and  only  a   small 
portion  is  the  recoil   energy   of   the  residual  nucleus   Y, 
one  can  consider  £~  to  be  the  kinetic   energy  of  b,    and: 


£f     =     &>Y        -         E^* 

where  £bY   is   given  by  £,  =  £bY   =£+   Qab   and  £  =   S*0 
is  the  kinetic   energy  in  the  entrance  channel    (    both  a 

and  X   are  in  their  ground  states).      This   represents  the 

maximum  value  of  C  which  can  occur  if   the  residual 

T 

8Blatt  and  Weiskopf :   344-345. 
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nucleus   Y  remains  in  the  ground  state.      Eo     ls  the  excita- 
tion  energy   of    Y   after   the  reaction  and   is   the   difference 
between  the   energy   of   Y   after   the   reaction  and  the 
energy   of   the  ground   state  of   Y.      So   long  as    E^     is 
less   than  £dY»    to    ever*y  lev©*   E/9*    of   the  residual   nucleus 

there   corresponds   an   energy  £    of   the   emerging  particle. 

r 

Should   the   energy  be  high   enough,    many   levels   of   the 
residual  nucleus   can  be   excited.      The   energy   distribution 
of   the  particle   b   becomes    continuous   when   the   energies 
Eg      are   closer   together   than  the   definition  of   the   energy 
in  the  incident   beam    or   closer   than  the   energy  resolution 
of   the  experimental  arrangement   for   detection  of   the 
reaction  products.       The   shape   of    the  distribution  function 
®b(£)    d£  f°r   the  particles   b    emitted  between  £  and  £  +d£ 
is   given  by: 

Gb(£)d£  =  lGc((3) 

where   the   sum   is    extended   over   the   channels/2  whose 
energy  lies   within  the   energy  interval   d£  .      The  number 
of   terms   in  this    sum   is   given  by   the  number   of   levels   of 
the  residual  nucleus   Y    with  an   excitation   energy   E/, 
between   E  and   E   -   d£   ,    where   E   =  cDY   ~  £  *      On©   can 
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call   this   number  W^(E)d£  and  WY(E)    the    "level   density". 
In  order   to   determine   the   relative  intensity   distribution 
of   the   emerging  particles   Ib(£)d£    ,    one  need   consider 
only  the   numerator   in  the   expression  for   the  probability 

GC^)J 

Ib(£)d£=    constant  ^c(^)wY(^bY   -£)d£ 
6~c(/9  )    =  6~c^)    IS   a   function  of   the   channel   energv  £  =  (fa  • 
The  factor      replaces   the  factor   k^2-   . 

Experimentally,    it   is   impossible   to   determine  the 
channels   by   which   a   particle  b   left   the  nucleus   or    even   to 
determine  its    energy.      One  is   therefore   interested  primar- 
ily in  the   cross-section  of   an    (a,b)    reaction  regardless    of 
the   specific    channels.      Should   the   energy   of   the   incident 
particle  be  high    enough   to   leave   the  residual   nucleus   in 
many   different   possible   excited   states,    one   can   express 
the  probability  in  terms    of   the  level   density  W.    Letting 
F^   =    Xkg,2-^^),    where   the   sum    extends    over   all   open 
channels B  leading   to   the   emission  of   b,    this    sum    can  be 
expressed   in  terms   of   an  integral: 

^bY 

Fb-  *V£bY)   -  ^     ^tfy)wY{£  bY  -£)d£ 

Qc(6&)    is    the   cross-section  for   the  formation   of   the 
10Blatt   and  Weiskopf:    367. 
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compound  nucleus   by  a    collision  with   the   energy  C -  between 
b   and  the   excited  nucleus   Y*,    the   latter  having  an   excita- 
tion  energy   ofX^Y   ~<Q"      These   cross-sections   depend 
only  upon  the   channel   energy/"  .      FD   is   a  function   of   the 

r 

maximum    energy  £^y    Whipk  the  particle  b    can  acquire  and 
leave   the  residual   nucleus   Y   in  its   ground   state.      If   the 
level   densities   Wy(0    are  known,    the  function  F^    can 
be   calculated   for   all   particle  types.       Thus   the   cross- 
section  for   an    (a,b)    reaction   can  be   expressed: 
^(a,b)    =  <TC  (a)    FK 

where   the   sum    extends    over  all   particle   types   c   that    may 
be   emitted  in  this   reaction. 

The  average   energy   E   of   a   system   is   a   monotonically 
increasing  function  of   the  nuclear   temperature  0.       The 
function   E(9)    has   a   vanishing   derivative/dE \  =o   for  0=o. 
If   one   expands    E(8)    near   9=o   in  powers   of   0,    it    must 
begin  at   least   with  a   quadratic   term.      This   leads   to   the 
following  assumptions: 

(1)  E   =    E(9)    permits   a   power   series    expansion 
around  0    =    o. 

(2)  the  leading   term    of   this    expansion  is   propor- 

Q- 

tional   to   0    . 


Blatt  and  Weiskopf:    369-370. 


(3)    the  higher   order   terms   of   9    can   be  neglected. 
One   can  now    write,    with  a  =    constant,    E   =   a  6,    from 
which: 

I       =  I    dE      =  2(aE)  +   constant i and  for 

/    e(E) 

the  level   density: 

W(E)    =    C    EXP(2lSE    ) 
This   is   a  statistical   formula   which   applies   to   a   region 
containing    many   energy   levels.      It   is   a   very  rough 
approximate   orientation  of   the   general  features   of   the 
densities   of   nuclear   levels.      These  level   densities    can  be 
measured   experimentally  by   comparing  the   energy   distribu- 
tion of   the   emerging   particles   in  a  nuclear   reaction  with 
the  results    obtained  from   the  theoretical  formula.    12 

d.      Secondary   Nuclear   Reactions 

In  most   nuclear   reactions   of   the  type  a   +   X   =Y   +   b, 
the   residual  nucleus   Y   is   left   in  an   excited   state.       The 
excess    energy  is    emitted  by   Y    many   times   in   the  form    of 
one   or    more  gamma-rays   until   such   time  that   the  nucleus 
reaches   the  ground   state.      In  some  instances   the   excita- 
tion of   Y    may  be  high    enough  for   Y   to    emit   another 
particle,    perhaps   a   neutron  or   a   proton.      One   can   say 


12Blatt  and  Weiskopf :    372 
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then  that  the  first  nuclear  reaction  is  followed  by  a 
"secondary  nuclear   reaction": 

Y    =   Z   +    c 
in  which  a  particle  c  is   emitted.      The  observed  reaction 
is   a  bombardment   of   X   by  a,    followed  by  an   emission  of 
b   and   c   leaving  a  residual   nucleus   Z,      One   can  represent 
this   process   as   an  X(a;bc)Z   reaction.      The   most   common 
reactions   of    this   type  are  the    (a;2n)    reactions.    *3 

Observed  X(a;bc)2   reactions    may  be   caused  by   two 
different   mechanisms: 

X    +   a   =   Y   +   b  Y   =   S   +    c  (A) 

X   +   a   =   Y»   +    c  Y'   =   Z   +   b  (B) 

If   the  time  between  the  emission  of  the  first  and  the 
second  particle  should  become  very   short    (of   the   order 
of   the  transit   time   of   nucleons   over  nuclear   dimensions), 
one   can  no   longer   refer   to   a   defined   excited   state   Y   or 
Yf   between  the  two    emissions.      The  reaction   should  then 
be   described: 

X    +   a-Z   +   b   +   c  (C) 

This   third  process    constitutes   a   simultaneous   breakup   of 
the   compound  system    C   into  a  nucleus   Z   and  two   particles 
b   and   c.      Process    (C)    does   not   correspond  to   any   of   the 
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"channels"   which  have  been  used   so   far.      Blatt  and 
Weiskopf   indicate  that   an  assumption  should  be   made   that 
the  reaction  X(a;bc)S   does   not   go   to   process    (C).    There 
is   some   doubt   as   to   whether   or   not   this   reaction  should 
be  included;    however,    in  the   computer   calculations    which 
follow,    this   assumption  is   not   made  and  the  process    (C) 
is  included.    ^+ 

The   cross-section  for  process    (A)    is    denoted  by 
£~  (a;bc)    and   the   cross-section  for   process    (B)    by 
5"  (a;c,b).       The   observed   cross-section   of   the    (ajbc) 
reaction   will   be   called  ^(a;bc).      It   is   the   sum   of    the 
cross-sections   for   processes    (A)    and    (B): 
^•(a;bc)    =  C(a;b,c)    +^(a;c,b) 
Usually   one  of   these  two   contributions   is   much   more 
important   than  the  other,    but   there   may  be  cases   in 
which  ^(a;b,  c)    and  c(a;c,b)    are  of  the  same  order   of 
magnitude.      The   distinction  between  processes    (A)    and 
(B)    should  not  be   made  if  b   and   c   are  the   same  kind   of 
particle    (e.g.,    (a;nn)    reactions).    *-> 

Considering  the   definition 

6^(a,b)    =  6~c(a)Gc(b) 

12j-Blatt  and  Weiskopf:    374-375. 
15Blatt  and  Weiskopf:    375. 
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the  cross-section   (a,b)    also  includes   those  reactions  in 
which   the  residual   nucleus   undergoes   a   secondary  reaction 
at  a  later  time.      One  can  therefore  distinguish  the  cross- 
section  £*( a, b)    which    constitutes   only  those  reactions 
where  the  residual  nucleus   does  not   emit  additional 
particles.     6~*(a,b)    is   what  is  understood   conventionally 
by   the   cross-section  of   an  X(a,b)Y   reaction.      The  follow- 
ing  relation  holds: 

6-(a,b)    =6~*(a,b)    +  £  «T(a;b,  c) 

c 

One  can  express  the  cross-section  for  the   (a;b,c) 
reaction   (process    (A))    as  follows: 

6^(a;b,c)    =6^  (o()  pGc^)GY/ff.   (c) 
0{  is   the   entrance   channel,    and  the  prime   on  the   sum 
indicates   that  it  is   extended  only  over  those   channels    B 
in   which  the  residual  nucleus   Y   is    excited   enough   to    emit 
a  particle   c.       Gc{fl)    is   the  relative  probability   (a  pure 
number)    that   the   compound  nucleus    C   will   decay  into 
channel  y^,    while  Gy*>(c)    is   the  relative  probability   that 
the  nucleus   Y   in  state^1    (corresponding  to    channel^) 
will   emit  a  particle  c.    10 

The  simplifying  assumption   can  be  made   that   the 
residual  nucleus   Y   always    emits   a   particle  if   its    excitation 


l6Blatt  and  Weiskopf :    376. 


13 


energy  is  high   enough.      Let  Smin(Y)    De  the  smallest 
separation   energy  of   a   particle  from   Y.      The  assumption 
is   then   expressed  by: 

Ic  GY/ff.(c)  ^  1  provided  that  £bY  -^  >    Smtri(Y). 

This   is   the   equivalent   of   neglecting   the   emission  of   gamma- 
rays   in   competition  with   the   emission  of   particles,    an 
assumption  which  is  approximately   correct   in   most   cases 
except  in  the  vicinity   of   the  threshold.      One   can  thus 
write: 

£-*(a,b)     ££"C(«V) 
where  the  double  prime   on  the  sum   indicates   that  it  is 
extended   only   over   those   channels      in  which  the  residual 
nucleus   Y   is    excited  below    smin(Y)»      Let  £sec  be  the 
maximum    energy  available  for   emission  of   secondary 
particles   from   the  first   residual   nucleus   Y   for   a   given 
excitation   energy  in  the   original   compound  nucleus.      This 
energy  is   also   the   energy  of   the   emerging  particles   b   in 
the  primary   X(a,b)Y   reaction  at   which  Y    can  just   begin 
emitting  particles    (See   Figure  1)  .  £  is   given  by: 

S  6C 

£sec   =   £bY   "   Smin(Y) 
One  then  obtains   the  following   expressions   for   the 

"observed"   cross-section  g"*( a, b)    of   the  X(a,b)Y   reaction: 


\h 
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Energy  Relationships  in  the  Reaction 
X+aZY*    +    bf 
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$-*(a,b)    =  ^(a,b)  (for  £         <0    ) 

c  usec 

•*bY 

6~*(a,b)    ^re(a,b)     ^sec (for<f>o) 

>bY 

7        ib(£")dr 


Where  1^,(6)    d£~     is   the   relative   probability  of   emitting 
particles   b   with   energies  £  between  f  and    /"+   d£  .      ' 

in 

(    Ib(<T)    d<f  =   constant(f^(/9)    WY(,fbY   -(f)    d£  )  .     9 

The   cross-section  for  a   secondary   nuclear   reaction 
with   emission  of   particles   c   following  a   primary    (a,n) 
reaction  can  be   written  with   the  same  approximation  as: 

L 


^(ajn,c)      ^r    £*(a,n) 


exp(    -  £-  )GY/9»   (c)d(f 

£bY 

exp(    -£)d£ 
o  « 


<f      =   £  -      S       (Y) 


In  this  formula  £_   is   the  excess    energy  over  the  thres- 

20 
hold  of  the    (ajn,n)    reaction.    SC(Y)    =   Sn(Y). 

Should   the   energy   of   the   incident   particle  be 

sufficiently  large,    the   residual  nucleus   after   emission  of 

the   second  particle    (nucleus   Z)    may  be   excited   enough   to 

17Blatt  and  Weiskopf  s    376. 

18Blatt  and  Weiskopf:    377. 
19Blatt   and  Weiskopf:    377. 
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emit  a   third  particle.      One   would  then  observe  tertiary 

21 
reactions   or   reactions   of    even  higher   order. 

3.         The   Statistical  Model 

The   model  portrayed  herein  is   a   statistical   representa- 
tion of   the  nucleus   and  has   been  adapted  for   calculations 
in  a   digital   computer.      This   model,    with   due  regard  to 
the  appropriate  assumptions,    is   general  and   can  be  used 
to    make   statistical   predictions   regarding  the   cross-section 
for   the   emission  of   particles   from  any  initial   nucleus.      In 
its   current   status,    the   computer  program    is   somewhat 
limited  in  that   the  following   nuclear   properties  have  been 
neglected: 

a.  The   coulomb   barrier 

b.  Spin 

c .  Parity 

d.  Angular   momentum 

e.  Tertiary  and  higher   order   reactions 
The   density   function   P(E)    is   the   same  for   the 

different  product   nuclei. 

Tertiary  reactions  were  considered,  but  have  not 
been  included  in  the  main  program .  There  is  included  in 
Appendix  II   a   computer   subroutine   which   can  be  used  to 
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compute  tertiary  reactions. 

The  coulomb  force   can  affect  the  neutron  processes 
indirectly.      Should  a   proton  process   be   made   smaller,    then 
the  neutron  process   will  be  made  larger.      In  the   calcula- 
tions  studied  herein,    the  coulomb   effects  are  neglected. 

In  setting  up  the  model  and  fitting  it  to  the  computer 
program,    Cu  ^   was   selected  as  the  experimental  nucleus 
with  which  to   compare  the  theoretical  results.      Cu"5 
was  selected  for  the  following  reasons: 

a.  Cu°5   is  a   light   nucleus   and  is  less   complicated 
than  others . 

b.  Cu  ^   has   a   low    coulomb  barrier. 

c.  Experimental   data   is   available  for   making 
desired  comparisons   with  theoretical   calculations. 

The   experimental  data  was  taken  from   graphs   of 
the  published  results  of   experiments   with  photo   neutron 
cross-sections   conducted  at  the   Ernest   O .    Lawrence 
Radiation  Laboratory,    Livermore,    California  in  August 
I963.    (See   Figures   2,3   and  k)    In  the   Livermore   experi- 
ments  the   single  neutron  cross-sections   actually   consist 
of  the  sum   of   the  reactions  s 

^n  +  Cnp      +  rn2P   + 

and  the  double  neutron  reactions: 

^~2n  +  r2np  +  ^2n2p  + 
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The  general  technique  here  is  to   compute  the  photo 
neutron   cross-section  for   the  X(Y,2n  )Y   process   and 
compare  it   with   the   experimental   data   for   the   same 
process.      This  is    an  application  of  the  relationship: 

J         exp(-|-)GY^(c)df 
^~(a;n,c)=  6"(a,n)    0 

fcbY 


f        exp(-|)df 


£    =  £         _   S    (Y) 

£~(a,n)    is   taken   from   the   experimental   data   for   Cu"^ 
and  is   used  throughout   in  all   computer   calculations.      For 
these  processes,   6"(a,n)    is  the  £{Ytn) ,    and  6(a;n,c)    is 
^(v}2n)  . 

The  range  of  gamma-ray  energy  for   which  calcula- 
tions  are   made   extends   from   0    to   30   MEV.      Within  this 
energy  range,    six  reactions  are  possible.      These  reactions 
are  portrayed  in   Figure   5  • 

In   order   to    make   computer   calculations   less   complicated, 
channel   widths   are   set  at   0.3   MEV,    and  all  integrals   are 
calculated  in  accordance   with  this  incremental  change  in 
energy.      Also,    the   threshold   energies   for   each   reaction 
have  been  adjusted   so   as   to   be   divisible  by  0.3   for   ease 
in  carrying  out   the   necessary  numerical   integration. 

The  probability   of   decay  to   any  particular    state  is 
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statistical,    and  is   computed  for  both   single-step  and  two- 
step   processes.      In   setting  up   the   model,    the  return  of 
the   compound  nucleus   to   its   ground   state  by   the   emission 
of   a   gamma-ray  is   neglected. 

From    Figure   5,    it   is   noted   that   there  are   two 
direct   processes   by   which   two   neutrons    can  be   emitted: 

a.  Cu  65  „  Cu63      +     2n 

b.  Cu^-^Ni62      +      p      +      2n 

For   simplicity,    the   Cu65    (Y;p,2n)Ni62   process   will  be 
omitted.       These   omitted   processes   contribute  to   the 
overall   cross-section   only  at   energies   above  their   thres- 
holds.     For   this   process,    the   threshold  is   taken  to   be 
23.7   MEV. 

Four   secondary   processes   result   in  the   emission  of 
two   neutrons  s 

a.  Cu^-».Cu6if      +      n — ►Cu63      +      2n 

b.  Cu6-^-*Ni624-      +      p — ►Ni62      +      p      +      2n 

c.  Cu6JL-*Ni63      +      p      +      n-^Ni62      +      p      +      2n 

d.  Cu6i^Cu63     +      2n-*Ni62      +      p      +      2n 
Two   tertiary   processes   result   in   the   emission  of 

two   neutrons  % 

a.  Cu^Ji- >.Ni6^      +      p_ ►Ni63      +      p      +      n->Ni62+   p   +   2n 

b.  Cu6-^->Cu6i+      +   n— ►Cu63      +2n-^Ni62      +      p      +      2n 
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In   each   case,    the  tertiary   processes   and   the   Cu°5 
(T;p,2n)Ni°^   processes   are  omitted   since   they   contribute 
to   the  overall   cross-section  only  at   energies   above   the 
threshold   of   23-7   MEV. 

Single  neutron   emitting  processes   are  as   follows: 
Three   direct 

a.  Cu6l_^cu64   +      n 

b.  Cu65 ►Ni63     +      p      +      n 

c.  Cu6^ ►Co62   +      2p      +      n 

The   (y;2p,n)    process   is   omitted. 

Three   secondary 

a.  Cu6^ ►Ni62*      +      p  ►Ni63   +      p      +      n 

b.  Cu6^ ►Ni6^      +      p ►Co62  +      2p  +      n 

c.  Cu65 ►Cu6^    +      n  ►Ni63   +      p      +      n 

The    (T;2p,n)    process   is   omitted.      All   tertiary   single 
neutron  processes   are  omitted. 

The  decay  of   the  compound  nucleus   for   each  process 

may   be   expressed  as   an  integral   of   the   density   function 

integrated   over  all  available   energy   channels: 

-E-Ethl 

P(E')dE' 
'0 

This   integral   represents   the   decay  of   Cu  ->   to    Nib^   +   p. 

Other   reactions   are  represented  accordingly. 
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Probabilities   for   direct   processes   are   calculated  as 
follows? 

a.  Two-neutron  processes 

(1)  Cu6l*Cu63      +      2n 

E-Ethl* 

p    (E)    m  y/    R(    E   -    Eth5    )    F(E')dE» 

E-Ethl         E-Eth2         E-Eth3         E-EthU         E-Eth$         E-Eth6 

P(E'dE»+  /p(EOdE,yP(EOdE'+/p(E')dE'+/P(El)dE«  +  ^(E')dE' 
o  'o  s0  'Q  JQ  SQ 

Where   P, (E)    is   the   probability   of   decay  of   the   compound 
nucleus   C*   to   Cu  ^   with   the   emission  of   two   neutrons. 
The  function  R(    E   -   Eth5    )    will   be   explained  in   subsequent 
discussion. 

(2)  Cu65 ►Ni62        +      p      +      2n 

E-Eth5 

P(E»)dE' 

P?(E)~  , l2_ 

E-Ethl    E-Eth2    E-Eth3    E-EthU    E-Eth5    E-Eth6 

P(E«  )dE'+/p(E»  )dE«+/p(E«  )dE«+/p(E»  )dE»  +/p(E«  )dE'  +  /p(E»  )dE» 
0       "0       'o  0       'o       '0 

b.  Single   neutron  processes 

(1)    Cu65 "-Cu64        +      n 

E-Eth2 

R(    E   - 

P3(E)=  «/q 

E-Ethl         E-Eth2         E-Eth3         E-EthU         E-Eth$         E-Eth6 

P(E»  )dE'+/ P(E«  )dE«yP(E'  )dE»+/p(E»  )dE'+/p(E»  )dE'+/p(E»  )dE» 

0  'o  '  Q  ' 0  -'O 
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(2)    Cu6-5 ►Ni63         +      p      +      n 


E-Eth3 


P,(E)=  ^) 


E-Ethl         E-Eth2         E-Eth3         E-EthU         E-Eth£         E-Eth6 
P(EOdE'+/p(E')dE'+/p(E')dE'-/ P(E«  )dE'+/p(E'  )dE'+/ P(E' )dE» 

0  -'O  -^o  -/o 

(3)    Cu  2 ►Co   '  +      2p      +      n      (omitted) 

The  probabilities   for   secondary  processes   are 
calculated   with  a   slightly   different   procedure.       These 
two-step  processes    consist   of  a   decay   of   the   compound 
nucleus   to   a   second   excited  nucleus.      This    second   excited 
nucleus   then   decays   by   emitting  a   particle.      Thus   the 
probability  for   the   complete   process   will  be  formed  by  the 
product   of   the  probabilities   of   two   independent   reactions. 
For   the  general    scheme   of   decay,    consider  a  hypothetical 
example.    (See   Figure   6) 

The  available   energy   of   the  initial   reaction  ranges 
between  0   and  30   MEV   and  is   represented  by   E.       Decay 
of  the   compound  nucleus   to    state  A   can   occur   if    E   is 
greater   than  the   threshold   energy   EtA.    Likewise,    if    E 
is   greater   than   EtB ,    the   compound  nucleus   may  also 
decay  to    state  B.      Should  this    condition  prevail,    then 
both   states  A  and   B    will   compete  for   the  reaction. 
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Probabilities   for   the   two-step   process    must   then   make 
provisions   for   the   consideration  of   this    competition.      Now, 
once  the  initial   reaction  has   taken  place  and  a   second 
compound  nucleus  has  been  formed,    this   second   compound 
nucleus   will   tend   to   decay  if   the  energy   is  high   enough. 
Thus,    if    E   -    EtA  is   greater   than   EtB    -    EtA,    this 
second   compound   nucleus   will   tend  to   decay.      The   decay 
of   this    second   compound  nucleus   by  the   emission  of   a 
particle   will   cause  a   decrease  in  the  probability   of   this 
second   state   existing   when  the   energy   exceeds   the   value 
of    E   -    EtB.      This    "leakage"    must   be   compensated  for 
when  the   probabilities   are   calculated. 

In   calculating   the  probabilities   for   two-step   reactions, 
the  reaction  which   occurs    second  physically   must   be 
computed  first   mathematically.      To    compensate   for   the 
"leakage"   factor,    the  function  R(    E   -   Eth)    is   defined 
as  follows; 

R(    E   -    Eth)    =    0      for      E   -    Eth  >  0 
R(    E   -    Eth)    -   1       for      E   -    Eth  <0 
The  particular   value   of   Eth   will   depend  upon  the   reaction. 
The  function  R(    E   -    Eth)    becomes   a   part   of   the   density 
function  for   the   second   reaction  and  is   integrated  along 
with   the   density  function. 
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The    "leakage"    must   also   be   considered  for   the   single 
neutron  and   single  proton  processes   for   energy   values 
above  the  two-neutron  threshold   since   the  two-neutron 
process   will    "rob"   the   single  neutron  process.       This 
"leakage"   is    compensated  for   in  the  appropriate   equations 
by  the   same  function   R(    E   -   Eth)  . 

Consider   the   two-step   reaction: 

X   — ►    Y        +        a ►Z        +      a      +      b 

Where   X   represents   the  initial   compound   nucleus,    Y 

represents   state  A,    and   2   represents   state   B.      Once 

state  B   has  been   formed,    this    state   will   exist   until   the 

energy   E2    exceeds   the  value   EtC-EtB ,    at   which   energy 

there  will  tend  to  be  a   "leakage"   of  particles  from   state 

B    to   state   C    (provided  of    course   that    state   C   is   one 

which  is   physically  possible  from    state   B).      This    "leakage" 

is   compensated  for   by   the  function  R(    E   -   Eth).       Thus 

the  integral   for   the   second   of   the  two   processes   would 

be: 

E   -    EtA   -    (  EtB    -   EtA) 

PAC(E)=    /  R(    El+EtA  -   EtB)    P(E')dE' 

0 

But,    E   -    EtA  -(EtB    -    EtA)    =    E   -    EtB,    and    El   +    EtA 

=    E.       Therefore: 

E   -    EtB 
Pac(E1)  =  /r(    E   -    EtB)    P(E')dE' 
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The  integral   for   the   complete  two-step   process,    combined 
with   the  appropriate   denominator,    would  appear   thusly: 


E   -   EtA 


POAc(E^^Er 
0 


E   -    EtB 


PAC(E1')    dEl'        dE' 


0 


D2 


Dl 

D  1  represents  the  denominator  for  the  decay  of  the 
compound  nucleus  to  state  A,  and  D2  represents  the 
denominator  for   the   decay   of   state  A   to    state  B. 

Using  the  above   procedure,    the   probabilities   for   the 
two-step   processes    can  be  written.      Consider  first   the 
secondary   reactions: 

a.      Ni62+*   +   p ►Ni63      +      p      +       n 

R(    E   -    Eth5)    P(E1»)    dEl' 
P0_  (  El )  = . _^___^__ 


2pn 


E   -   Eth3  E   -   Eth5  E   -    Eth6 

P(El')dEl'   +  /p(El')dEl'   -f  /  F{  El' )  dEl' 


,         ^  64*  ^  .T.63 

b.      Cu  ^     +  n 


P2np<E2>  = 


-►  Ni  +      n      +      p 

R(    E   -    Eth5)    P(E2')    dE2' 

E   -    Eth3  E   -    Eth5  E   -    Eth6 

P (  E2« )  dE2'+  /  P(  E2' )  dE2'  +  /p (  E2» )  dE2' 
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Now,    the  probabilities   for   the   two-step   single-neutron 
processes    can  be   written: 

a.  Cu65^ ►  Ni6^*      +   p— ►Ni63      +   p  +   n 

E   -    Ethl      E   -    Eth3 

^-)  -  m 

P        (  E )  - ^_0 SO 

J3&S1         E-Eth2         E-Eth3         E-EthU         E-Eth5         E-Eth6 

P(E»  )dE« ■»/  p(E»  )dE'+/P(E'  )dE»+/ P(E«  )dE'+/P(E»  )dE'+/p(E'  )dE' 
0  ^0  ^0  "  0  ^0  '0 

b.  Cu^* ►cu6^      +  n — ►Ni63      +  n  +  p 

E   -    Eth2      E   -    Eth3 
P(E') 


P0       (E2)    dE2    dE» 
2np 


P2np<E>  = 


E-Ethl         E-Eth2         E-Eth3         E-EthU         E-Eth$         E-Eth6 
P(E«  )dE»+/p(E«  )dE»+/p(E«  )dEfy  P(E«  )dE,-/p(E»  )dE«+/ P(E»  )<JE» 


The  probability  for   the   two-step   two-neutron 
process   can  be   obtained   similarly.      Again  consider   first 
the   secondary  reaction: 

Cu6^*   +   n-^Cu63      +   2n 


P2nn(E2)=s  R    (     E   -    Eth5)    P(E2')dE2' 

E-Eth3  E-Eth4  E-Eth5 


E-Eth6 
E2'+  fF(E2')dE, 
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The  probability  for   the   complete   two-step   process 
becomes : 

Cu65*-^  cu6j+*      +      n— *-Cu63      +      2n 

E   -   Eth2  E   -   Eth4 


P2nn<E)  = 


0 


0 


E-Ethl    E-Eth2    E-Eth3    E-EthU    fi-Eth£    E-Eth6 
P(E«  )dE»+/P(E«  )dE«+/p(E»  )dE'+/ P(E«  )dE'+/p(E'  )dE«+/>(E'  )dE» 
0       ^0      ^  0      "0        0      '  0 

The   total   probability  for   each   process   is   the   sum   of 
the  probabilities   for   the  individual   reactions   since   these 
are   mutually   exclusive   events.       The  total   probability  for 
the  two-neutron   processes    can  be   written: 

Pnn<E)    -   P2nn<E)    +    ^E>    +    P2<E> 
The   total   probability  for   all   single-neutron  processes   is: 

Pn(E)    =    P2pn(E)    +    P2np(E)    +    P3(E)    +    P4(E) 
Statistically,    the   ratio   of   the   cross-sections   for 
the  two-neutron  processes   to   the   cross-sections   for   the 
single-neutron   processes    should   equal   the  ratio   of    corres- 
ponding probabilities: 


£  (Y",2n) 
6"  (V,n) 


P(V,2n) 
P(V,n) 


and 


6~CT,2n)    =^(V,n) 


P(V,2n) 
P(V,n) 
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where   P(V,2n)    =   Pnn(E)    and   P(T,n)    =   Pn(E)'       This 
relationship   is   the  basis   on  which   the   computer   program 

is   written. 

2f.        Computational   Parameters 

Blatt  and  Weiskopf   present   the  function 

P(E)    =   C    exp(27AE    ) 
as  a  first   guess   in  the  determination  of   the   density- 
function  describing   the   statistical   decay   of   the   compound 
nucleus.      In  this   function,    E   represents   the   energy  of 
the  incident  particle   (gamma-ray);    C  is  a   constant   of 
proportionality;    and  A  is   a   constant   based   on  the  nuclear 
temperature  0 . 

From   this  first  guess,    several  variations  have  been 
devised  in  an  attempt   to   obtain  a   density   function  which 
will   closely  approximate  the  available  experimental  data 
for  the    (V,2n)    process.      In  each  of   these  density  functions 
the  following  constants  are  used: 

A  =  Nuclear  temperature  constant;    the  units   of 
A  are    (MEV)'1 

C   =   Constant   of   proportionality 

22Fultz,  S. C, Bramblett, R. L.  ,  Caldwell,  vJ.T., 
and  Harvey,  R.R.,  Photo-Neutron  Cross-sections  for 
Natural  Cu,  Cu^3,  and  Cu&5.  Physical  Review,  v.  133, 
no.    5B,9  March  I964:    B1149   -  B115^. 

23Fultz,    Bramblett,    Caldwell  and   Harvey: B1149-B115i+. 

22*Fultz,    Bramblett,    Caldwell  and   Harvey :B1149-B1151+. 
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n  ■   RN   =   Exponent   of   the    (AE)    term 
Altogether,    fourteen   density  functions  have  been 
formulated  for  this   program    (See  Appendix  III);    however, 
only  six  have  been  investigated  in  any  detail.      These  six 
ares 

a.  P(E)    =   C(S   -    E)    EXP(A(En   -   Zn)  )  J    2   =    E-Eth 

b.  P(E)    =   C    EXF(    A(En)) 

c.  F(E)    -   CZ    EXP(    A(    En)) 

d.  P(E)    -    CE    EXP(    A(    En)) 

e.  P(E)    -   CE/EXP(A    (    En)  ) 

f.  P(E)    =   C(E)3    EXP(A    (En)) 

Each  of  these  density  functions    has     been  tested 
•with  various  values   of  A.      In  all   cases,    the  parameter 
n  =   RN   was   set    equal  to   1.0.      The  values   of  A  ranged 
between  0.0001  and   10.0      (MEV)"1. 

For   each  of   t"he  values   of  A,    the  cross-sections  for 
the  theoretical    CK,2n)    process  were  plotted  versus   the 
energy   of   the  incident   gamma-ray.      On  the   same  graph  in 
each   case  were  plotted  the  experimental   cross-sections  for 
the  corresponding  values   of   gamma-ray  energies.      The 
procedure   used  is   to    compare  the   theoretical   cross-sections 
with   the   experimental  values.      The   objective  here   is   to 
find  a   density  function  P(E)    that   will   give  a   close 
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approximation  to   the   experimentally   determined   cross- 
sections. 
5.         Results 

a.  P(E)    =   C(3   -   E)    EXP(    A(    En-Sn)  }    Z=E   -Eth. 
See  Figures    7   through   14 .      The  idea   for   this   function 
was   obtained  from    Ericson's   article  on  The   Statistical 
Model  and  Nuclear   Level   Densities.      At  values   of  A 
ranging  from   0.0001   to   1.0    (MEV)     ,    it   shows   a   shape 
that  in  general   conforms  to   the  shape  of   the   experimental 
curve;    however,    for   larger   values   of  A  the   curve   develops 
a   spike  near  the  lower  energy  values.      Thus,    this  function 
is  not   considered  to  be  appropriate. 

b.  P(E)    =   C    EXP(A(En)).      See   Figures   15 
through  23  •      This   function  is  a   modification  of   the   one 
suggested  by  Blatt   and  Weiskopf  as   a   first   guess.      This 
function  was  tested  for  values   of  A  ranging  from   0.0005 
to   10.0    (MEV)"  .      This   function  reached  its   maximum 
value  at  an  A  value  of  about  0.1    (MEV)""1  and  held  this 
value  for  all  decreasing  values   of  A.      As   the  value  of  A 
increased  from   0.1    (MEV)      »    the  function  decreased  in 
value.      The  values   became  so   small  at  an  A  value  of  about 
2.0    (MEV)"1   that   the   curve   did  not    even   show   up.      In 

all  instances,    a   small   spike  appeared  on  the  lower   energy 
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end  of  the  curve |    however,    it  was  not  as  pronounced  as 
the  spike  obtained  with  the  function  in  a.    above.      This 
function  is   not    considered   to   be  descriptive   of   the 
experimental  cross-section. 

c.  F(E)    -   C2    EXF(A   En)  .      See   Figures   24 
through  33  •      For   this   function,    A  ranged  in  value  from 
0.0001   to   10.0    (MEV)"1.      Throughout   the   entire   range  of 
A  values,    the  shape  of  the  curve  was  generally  in   conform- 
ity  with  the  experimental   curve  except  that  the  ordinates 
were  at  all  times   much  too   small.      This  function  reached 
its   maximum   value  at  an  A  value  of  about   0.5    (MEV)"^. 
For  increasing  values,    the  curve  became  unobservable  at  an 
A  of   about   1.5    (MEV)"^.       This   function  gives   a  good 
shape,   but  the  values   of  the  cross-section  are  too   small. 

d.  F(E)    =   CE    EXF(    A(    En)  )  .      See   Figures  3k 
through  49.      Values   of  A  ranged  from   0.0001  to  10.0 
(MEV)     .      The  shape  of  the   curve  was  good  throughout, 
but  its   maximum  value  never   exceeded  about  one— fourth 
of   the  experimentally  determined  values.      It  reached  its 
maximum  value  at  an  A  of  about  0.01    (MEV)-1;    it  became 
too   small  to   observe  at  an  A  of  about  1.5    (MEV)     • 
T^his  function  does   not  have  a  high    enough    maximum   value. 
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e.  P(E)    =   CE/EXP(    A(En)).      See   Figures    50 
through   56.      The   values   of  A  ranged  from   0.0001   to 
10.0    (MEV)"1   for   this   function.      Throughout   this   range 
of  values,    the  curve  was  visible  and  had  significant 
magnitude.      However,    there  developed  a   spike  on  the  low 
energy   end  which  became  quite  pronounced  as   the  values 
of  A  increased.      This  is   not   considered  to  be  an  approp- 
riate  density  function. 

f.  F(E)    =   C(E)3    EXP(    A(    En)).      See  Figures   57 
through  64.      This   function  was   devised  in  an  attempt  to 
increase  the  magnitude  of   the  ordinates  for  the  cross- 
section   curve.      From   the  figures,    it   can  be   seen  that 
instead  of  increasing  the   magnitudes   of  the  ordinates,    it 
actually  decreased  them.      For  this  function,    values   of  A 
ranged  from   0.0001   to  3.0    (MEV)"1.      Throughout    most   of 
this  interval,    the  curve  held  a  good  shape.      The  magnitudes 
became  so   small  that  as   the  values   of  A  increased  through 
1.0    (MEV)-1   the   curve   was   no   longer   visible.      Although 

the   energy   coefficient    (E)3    does  not  adequately   describe 
the  process,    it   does    suggest   the  possibility  that   perhaps 
some   other   power   of    (E)    may   give  the   curve  a  better 
shape . 
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6.        Conclusions 

None  of   the  density  functions   tested  thus  far  give 
an  adequate  approximation  to  the  experimentally  obtained 
data.      Two   of   the  tested  functions   produce  a  good  shape, 
but  in  all   cases  the  maximum  values  do   not  approach 
the  experimental  values.      The  lack  of    conformity   so  far 
may  be  due  to   several   causes: 

a.  The  proper  density  function  has   not  been  found. 

b.  Coulomb  and  other   effects  have  been  ignored. 

c.  Tertiary  reactions  have  been   omitted. 

The  results   indicate  that   there   may  be   some  power 
function  of    (E)    which  acts  as   a  variable   coefficient   of 
the   exponential   term.       The  density  function  with   the    (E)3 
coefficient  is   a  first  attempt   to   find   such  a  term. 

The   results    obtained,  although   not   in   close  agreement 
with   experimental  data,    are   encouragng  and  indicate  that 
this    model   with   appropriate    modifications   and  development 
can  be   used  for    making  useful   statistical   approximations 
for  niclear   reactions. 
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APPENDIX    I 
CROSS-SECTION   PROGRAM 
This   program    was   written  for  the  purpose   of 
calculating   cross-sections  in  the  decay  of   the  compound 
nucleus.      Although  written  for  the  X(Y,2n)Y   reaction 
with   Cu^  as   the  target  nucleus,    this  program    can  be 
used  for  any  nucleus  in  computing  an  approximation  to   the 
cross-sections   for   the   emission  of  particles.      Appropriate 
assumptions   must  be  made  with  regard  to  the  coulomb 
barrier,    spin,    angular   momentum,    etc. 

In  general,    the  program   consists   of  a   series   of 
calculations  in  which  three  subroutines   perform   most   of 
the  mathematical  operations.      These  subroutines   calculate 
values   of   single  and  double  integrals.      A  subroutine  for 
the  computation  of   triple  integrals  is  included,    but   was 
not  used  in  the  main  program. 

Subroutine  PROD   calculates  the  values   of   single 
integrals.      The  call  statement  is: 

CALL    PROD(    EL,    PIN,    PNO,    NWB ) 
The  variables   in  the  call   statement  have  the  following 
meaning: 

EL    is   the   threshold  for   the  reaction  being   calculated, 
PIN  is   a   constant  used  in  calculating  double  and 
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triple  integrals.    For   the   single  integral,    PIN   equals   1.0. 

NWB    is    the  index   corresponding  to   the   energy   range 
over   -which  the   calculation  is    made.      It   is   determined  by 
dividing  the   energy  in  MEV   by   0.3-    (e.g.,    NWB    =    E/.3). 

It  is  appropriate  to   note  here  that  all  threshold 
energies  have  been  adjusted   so   as   to   be  divisible  by   0.3 
in  order   to   simplify  the   computer   operations. 

Subroutine   GRASP   calculates   the  values   of   double 

integrals.       The   call   statement   is: 

CALL    GRASP    (  ELI,  EL2,  EL3,  EL4,  EL 5,  EL6, 

PN20,POCN,PNLN,EFU,NPT,NP1,NP2,NP3,NP4,NP5, 

NP6). 

In  order   to   observe   more   clearly   the   meaning  of    each 

variable,    consider  a   two    step   reaction: 

X(Y,n)Y*  ;      Y* — ►S      +      n 

For  the  total  reaction  of  X   +Y — ►  Z  +  2n  over  the 

energy   range   of   0    to   30   MEV,    six  reactions   are  possible. 

Let   the  probability  for   the  first   reaction  be: 

E-Ethl 

P(E')P2ndE' 

P,    (E)=  '  0 

ln  E-Ethl X=m2 X=Em E-EthU E^^H3 E-Eth6 

P(EO<ffit+/p(EOdE'+yP(E0^,+/P(EOdE'yP(E0^t+/p(E,)<iEl 

o  'o  ^o  y o  '  o  Jo 
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Now   PNLN   =   PNL1   corresponding  to   the   Ethl  for   the 
threshold   one  reaction  in  the   numerator.       ELI   through    EL6 
are  the  threshold   energies   for    each  reaction  that   can  take 
place   over   the  prescribed   energy   range.       EPL    is   the 
threshold   energy  for   the  initial   reaction;    in  this    example 
EEL   =    ELI. 

Next,    consider   the   second   reaction: 

E   -   Eth3 
/r(    E   -    Eth5)    P(E3')    dE3' 

P2n<E>    "  2 

E   -    Eth3  E   -    Ethi+  E   -    Eth6 

P (  E3' )  dE3'+     /  P  (  E3' )  dE3'+  /  P  (  E31 )  dE3' 
0  -/o  ^0 

In  the   call   statement,    POCN   =   POC3    corresponding  to 
the   Eth3   for   the   threshold   three  reaction  in  the  numerator, 
EPU   is   the  threshold  for  the   second   reaction  and  in  this 
example    EPU=   EL3.      The  variables   NP1  through   NP6  are 
indices   for   the  computation  of   the  integrals   in  the  denom- 
inator.     If   a   particular   integral   is   not  to   appear   in   the 
denominator,    then  the   corresponding  index   will  be   set 
equal  to  2;    if   the  integral  is  to  appear  in  the  denominator, 
the  index   will  be  1.      For    example,    in  the  probability  for 
the   second  reaction  the  integral   for    Ethl   does   not   appear 
in  the   denominator;    therefore,    NP1   will   equal  2   in  the  call 

statement.      For   the  above  probability,    the  sequence   NP1 
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through   NP6    will  appear   in  the   call   statement   as   2,2,1,1,2,1. 

The  value  for   NPT  is  the   "leakage"   factor.      In  this 
example   NPT    =    Eth5/0.3. 

PN20    is   the   output   of   this   subroutine  and  it 
represents   the   probability   of   the  formation  of    cross- 
sections   by  a   two-step   process.      PN20    may  have  any 
name  in  the   call   statement. 

The  triple  integral   subroutine  follows   the   same 
general  procedure. 
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APPENDIX  II 
PROGRAM  FLOW  CHART 


START 


I 


READ  SEA 

FOR  ALL  n 
PROCESSES 


I 


READ  SEC 

FOR  ALL  2n 
PROCESSES 


I 


E(I)  =0.0 


I 


E(J-hl) 

=  E(J)  +  0.3 

NO 


"  YES 
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COMPUTE 
PNC1 
SUBROUTINE 
PROD 


I 


COMPUTE 
PNC2 
SUBROUTINE 
PROD 


COMPUTE 

PNC3 

SUBROUTINE 

PROD 


COMPUTE 
PNC*f 
SUBROUTINE 
PROD 


COMPUTE 

PNC5 

SUBROUTINE 

PROD 


COMPUTE 
PNC6 
SUBROUTINE 
PROD 
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DENOM(J)  =  PNC1  -I-  PNC2  +  PNC3 
+  PNC*+  +■  PWC5  +  PNC6 


J<?7 


J>57 


PNCX      = 
PNC2(I) 


PN6(J)    = 
PNCX/ 

DMOM(J) 


i 


PNCX      ~ 
PNC2(J) 
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YES 


K  =  J 


J>79 


.T<^7Q 


PNCY      = 
PNC3(K) 


PN7(J)  = 

PNCY/ 
DENOM(J) 


1 


PNCY      = 
PNC3CJ) 
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-   J>79 


J<79 


PN>+(J)    = 

PNCZ/ 
DENOM(J) 


] 


PNCZ       - 
PNC*+(J) 
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COMPUTE 

pro 

SUBROUTINE 
GRASP 


COMPUTE 
PNP 
SUBROUTINE 

GRASP 


I 


COMPUTE 
PNU 
SUBROUTINE 
GRASP 


I 


PNJ(J)  =  PN^(J)  -•-  PNU(J) 


I 


PNK(J)  =  PN6(J)  +  PN7(J)  + 
PNR(J)  +■  PNP(J) 


I 


CRN(J)  - 


SEA(J)  PNJ(J)/ 
PNK(J) 
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PLOT 

CRN (J) 

SUBROUTINE 

DRAW 


I 


PLOT 
SEC  (J) 

SUBROUTINE 
DRAW 


END 
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SUBROUTINE     PROD 


START 


I 


K  =  l 


m  -  *3 

Z=E(K)-EL  j 

J 


I 


W^LO.(Z/3.) 
NW=Wf  U 


Y=0„0 


I 


P(L)=F(Y,Z) 


I 
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no 


sum  no»o 


NY-W 
L  =  2 


SUM  I  =  SUM  I  +  P(L) 


NO 


111 


PNO(K)  "  H/2.(  P(l)  +  2.(  SUM  I  ) 

+-  P(NW)) 


<3 


NO 


KzK-hl 


RETURN  TO 
MAIN 
PROGRAM 


C     END   J 
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SUBROUTINE  PROBE     II 


START 


NWB  -1 


YES       /  IS 

E(NWB)<0; 
? 


NO 


H  3#3 
ZrE(NWB)-EL 


I 


W=10*(Z/3O 

NW  =  W  +  1 


Y=0.0 
L=  1 


P(L)ZF(Y,Z) 
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SUM  I  -0.0 


I 


L  -2 


I 


SUM  1=  SUM  I  4  P(L) 


NO 


114 


RETURN 
TO 
MAIN 
PROGRAM 


i 


END 


3 
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SUBROUTINE  GRASP 


START 


J  =  1 


YES 


NWB  ~  10.    (E(J)    -  EPL)/3< 

+■  1.0 


® 


EP  ~  E(J)   -  EFU 


,»   NO 


PIN  -  1. 


YES 

— ►- 


<s> 
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YES 

/IS  \ 
<T    EL<0    > 

?   NO 

POCl  =0.0 

NO 

/computer 
\poci  / 

f   YES 

\ 

> 

COMPUTE 
POCl 
SUBROUTINE 
PROBE  II 

t 

ELZEL2-EPL 

YES 


P0C2  ~   0.0 


COMPUTE 
P0C2 

SUBROUTINE  i 
PROBE  II 
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YES 

/  IS  \ 
<    EL  <0    N 

[NO 

P0C3  -  o.o 

NO^, 

>/COMPUTE\ 
\  P0C3  / 

\y/ 

f  YES 

i 

i 

COMPUTE 

P0C3 

SUBROUTINE 
PROBE  II 

t 

ELZEL^—EPL 

COMPUTE 
POO1*- 
SUBROUTINE 
PROBE  II 
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YES 

/    IS    \ 

<.    EL  <0     > 

T  NO 

poc5=o»o 

NO 

/COMPUTE rs 

\poc5  / 

J    YES 

i 

1 

COMPUTE 

POC5 

SUBROUTINE 
PROBE  II 

♦ 

EL~  EL6~EFL 

YES 


POC6rOeO 


COMPUTE 
POC6 
SUBROUTINE 
PROBE  II 


11 


J  >NPT 

-*- 


J<NPT 


POCPrPOCN 


I 


POCN  ZPOGP 


PNA  r  POCN/(  POC1  +■  P0C2  +■  POC3 
+  FOGh  -h  POC5  +  POC6  ) 


PNA  =  0.0 
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COMPUTE 
PNL1 
SUBROUTINE 
PROBE  II 


COMPUTE 
PNL2 
SUBROUTINE] 
PROBE  II 


COMPUTE 

PUL3 
SUBROUTINE 
PROBE  II 


COMPUTE 
PNL*f 
SUBROUTINE 
PROBE  II 


COMPUTE 
PNL5 
SUBROUTINE 
PROBE  II 


COMPUTE 
PNL6 
SUBROUTINE 
PROBE  II 
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10 


FN20(J)  ~  PNLN/(  PNL1  4-  PNL2  + 
PNL3  +  PNlA  -f  PNL5  +■  PNL6  ) 


PN20(J)=0 


RETURN 
TO 
MAIN 
PROGRAM 


END 
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APPENDIX    IIA 


SUBROUTINE  GROPE 


YES 


NWC  r  10o(  E(J)  -  EPL  )/3o  +  1. 


YES 


NWBr  IO»(  E(J)  -  EPL1  )/3.  +■  1.0 


0 


PIN  ~  1.0 


I 


EL  -  KL1  -  EPL1 


YES 


PN3A1-0.0 


COMPUTE 
PR3A1 
SUBROUTINE 
PROBE  II 


I 


EL-  EL2  -  EPL1 
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COMPUTE 
PN3A2 
SUBROUTINE 
PROBE  II 


I 


ELz  EL3  -  EPL1 


125 


COMPUTE 

PN3A3 
SUBROUTINE 
PROBE  II 


I 


EL  Z  EL*+  ~  EPL1 


YES 


PN3A*f=:  0.0 


COMPUTE 

? 

YES 


COMPUTE 
PN3A^ 
SUBROUTINE 
PROBE  II 


I 


EL-  EL5  -  EPL1 


I 
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YES 


PK3A5zo*o 


COMPUTE 
PN3A5 
SUBROUTINE 
PROBE  II 


I 


ELZ  EL5  -  EPLL 


YES 


PN3A6~0.0 


COMPUTE 
PN3A6, 

YES 


COMPUTE 
PN3A6 
SUBROUTINE 
PROBE  II 
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12 


J  <NPQ 


FET3A  =  PN3M/(   PN3A1  4-  PN3A2  4-  PN3A3  + 
PN3A*f  +■  PN3A5  +  PN3A6 


PN3A  =0.0 


PINU  PN3A 


I 


ELZ  ELI  -  EPL 


t 
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YES 


poci-  CUO 


I 


POBl  =  0-0 


COMPUTE 
POCI 
SUBROUTINE 
PROBE  II 


COMPUTE 
POB1 
SUBROUTINE 
PROBE  II 


I 


EL  =  EL2  -  EPL 
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YES« 


P0C2  =  0.0 


P0B2  =  0.0 


COMPUTE 
P0C2 
SUBROUTINE 
PROBE  II 


COMPUTE 
P0B2 
SUBROUTINE 
PROBE  II 


EL  Z  EL3  -  EPL 
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YES 


P0C3  =  0.0 


I 


P0B3  =  0.0 


COMPUTE 
P0C3 
SUBROUTINE 
PROBE  II 


COMPUTE 
P0B3 
SUBROUTINE 
PROBE  II 


EL  —  Elk   -  EPL 
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YES 


POC&  -  0.0 


I 


POB*f  =  0.0 


COMPUTE 
POC^ 
SUBROUTINE 
PROBE  II 


COMPUTE 
POB*f 
SUBROUTINES 
PROBE  II 


I 


EL  ~   EL5  —  EPL 


132 


COMPUTE 

poc5 

SUBROUTINE 
PROBE  II 


COMPUTE 

pob5 

SUBROUTINE 
PROBE  II 


EL  —   EL6  -  EPL 


133 


YES 


P0C6Z  0,0 


P0B6  Z  0.0 


COMPUTE 
P0C6 
SUBROUTINE 
PROBE  II 


COMPUTE 
P0B6 
SUBROUTINE 
PROBE  II 
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J>NPR 


J<NPR 


YES 


POCP-  POCN 


I 


POCN=POCP 


PN1A1  -  POCN/(  POB1  -h  P0B2  ■+-  POB3  + 
POB^  H-  POB5  +  POB6  ) 


PINHPNIAI 


t 
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COMPUTE 
PNL1 
SUBROUTINE 
PROBE  II 


COMPUTE 
PNL2 
SUBROUTINE 
PROBE  II 


COMPUTE 
PNL3 
SUBROUTINE 
PROBE  II 


I 


COMPUTE 
PNlA 
SUBROUTINE 
PROBE  II 


COMPUTE 
PNL5 
SUBROUTINE 
PROBE  II 


COMPUTE 
PNL6 

SUBROUTINE 
PROBE  II 
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FN1(J)  ~  PNLN/(  PNL1  ■+-  PNL2  H- PNL3 
+  PNlAf  -h  PNL5  +  PNL6  ) 


RETURN 
TO 
MAIN 
PROGRAM 


C 


END 
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APPENDIX    IIB 

PROGRAM  LISTING 
PROGRAM  PRODUCE 


DIMENSION  SEA( 101)  , SEC (101) ,  Y ( 102 ) , I T I TLE ( 12 ) , 

1  El  102)  »PNU( 101) ,PN4( 101) tPN5( 101)  ,PNJ( 101) » 

2  PNCK 101) , PNC2(101>  ,PNC3( 101) ,PNC4( 101 )  , 

3  PNC5( 101) ,PNC6(101) ,PN7( 101 ) ,PN6( 101) tCRN< 101) , 

4  PNR(lOl) ,PNK( 101) »PNP(101 ) ,DENOM( 101) 
COMMON  E»POCl,POC2,POC3»POC4,POC5,POC6, 

1  PNL1,PNL2»PNL3,PNL4,PNL5»PNL6,RN,'C,A> 

2  DENOM 
RN  =  1.0 
C  =  1.0 
A  =  1.0 
READ  99  SEA 
READ  99  SEC 

99  FORMAT  (  20F4.0) 
E(l)  =  0.0 
DO  50  J  =  1,101 

50  E(J  +  1  )  =  E(J)  +  .3 
NWB  =  101 

PIN  =  1.0 
EL  =  7.2 

CALL  PROD  (  EL,PIN,PNC1,   NWB  ) 
EL  =  9.9 

CALL  PROD  (  EL,PIN»PNC2»   NWB  ) 
EL  =  17.1 

CALL  PROD  (  EL,PIN,PNC3,   NWB  ) 
EL  =  17.7 

CALL  PROD  (  EL,PIN,PNC4,   NWB  ) 
EL  =  23.7 

CALL  PROD  (  EL,PIN,PNC5,   NWB  ) 
EL  =  26.7 

CALL  PROD  (  EL,PIN»PNC6,   NWB  ) 
DO  6161  J  =  1,101 
PNC22(J)  =  PNC2(J) 
PNC33(J)  =  PNC3(J) 
PNC44(J)  =  PNC4U) 
6161  CONTINUE 

DO  990  J  =  1,101 

DENOM(J)  =  PNCl(J)  +  PNC2U)  +  PNC3(J)  + 
1  PNC4(J)  +  PNC5(J)  +  PNC6(J) 
IF(  J  -  57  )52,49,51 
49  I  =  J 

51  PNCX  =  PNC22C I ) 
GO  TO  57 

52  PNCX  =  PNC2(J) 
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57  PN6(J)  =  PNCX/DENOM( J) 
IF(  J  -  79  )59,53,54 

53  K  =  J 

54  PNCY  =  PNC33(K) 
GQ  TO  67 

59  PNCY  =  PNC3(J) 
67  PN7(J)  =  PNCY/DENOM( J) 
I  F.I  J  -  79  )103,101.102 

101  M=  J 

102  PNCZ  =  PNC44IM) 
GO  TO  77   • 

103  PNCZ  =  PNC4(J) 

77  PN4(J)  =  PNCZ/DENOM( J) 
990  CONTINUE 

CALL  GRASP  (  7.2,9.9,17.1,17.7,23.7,26.7, 

1  PNP,POC3,PNL2, 17. 1,9.9.79, 

2  2,2,1,2,1,1  ) 

CALL  GRASP  (  7.2,9.9,17.1,17.7,23.7,26.7, 

1  PNR,P0C3,PNL1, 17. 1,7.2,79, 

2  2,2,1,2,1,1  ) 

CALL  GRASP  (  7.2,9.9,17.1,17.7,23.7,26.7, 

1  PNU,POC4,PNL2,17.7,9.9,79, 

2  2,2,1,1,1,1  ) 

DO  1000  J  =  1,101 
PNJ(J)  =  PN4(J)  +  PNU(J)   ^ 

PNK(J)  =  PN6(J)  +  PN7(J)  +  PNR(J)  +  PNP ( J ) 
CRN(J)  =  SEA(J)  *(PNJ(J))/  PNK(J) 
1000  CONTINUE 

NUMPTS  =  .101 
MODCURV  =  1 
Y(l)  =  100.0 
DO  41  J  =  2.101 

41  Y(J)  =  CRN(J) 
ITYPE  =  0 
LABEL  =  4H 
ITITLE(I)  =  8HLOWRY 
ITITLE(2)  =  8HA  =  1.0 
ITITLEO)  =  8HRN=  1. 
ITITLEI4)  =  8HE=ENERGY 
ITITLE(5)  =  8H1  UNIT 
ITITLE(6)  =  8H=  1  MEV 
ITITLE(7)  =  8HY=CROSS 
ITITLE(8)  =8HSECTION 
ITITLE(9)  =  8H1  UNIT 
ITITLE(IO)  =8H=  10  MB 
ITITLE( 11)  =  8H   F  -  1 
ITITLE(12)  =  8H  CU  65 

CALL  DRAW  (  NUMPTS ,E ♦ Y  , MODCURV , I  TYPE , LABEL , I T I TLE , 

1  EXSCALE,  YSCALE,  IXUP,  IRIGHT,  MODEXAX,  MODEYAX,  IWIDE, 

2  IHIGH,  IGRID,  LAST  ) 
MODCURV  =  3 

DO  42  J  =  '2,101 

42  Y(J)  =  SEC(J) 

CALL  DRAW  (  NUMPTS ,E , Y , MODCURV , I  TYPE , LABEL , I T I TLE , 
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1  EXSCALE,  YSCALE*  IXUP»  IRIGHT,  MODEXAXt  MODEYAX,  IWIDE 

2  IHIGH,  IGRID,  LAST  ) 

END 

SUBROUTINE  GRASP  (  EL1,EL2,EL3, 

1  EL4,EL5»EL6»PN20,POCN»PNLN»EPU,EPL,NPT, 

2  NP1,NP2»NP3»NP4,NP5»NP6  ) 
DIMENSION  E( 102) ,PN20( 101 ) ,DENOM(101  ) 
COMMON  E»POCl»POC2»POC3»POC4»POC5»POC6, 

1  PNL1,PNL2»PNL3,PNL4,PNL5»PNL6»RN,C»A» 

2  DENOM 

DO  889  J  =  1»101 
IF(EU)  -  EPL  )10,3,3 

3  CONTINUE 

NWB  =  10.*((E(J)  -  EPD/3.)  +  1.0 
EP  =  E(J)  -  EPU 
IF(EP) 11,4,4 

4  CONTINUE 
PIN  =  1.0 

1  EL  =  ELI  -  EPL 
IF(EL)22.12,12 

12  GO  TO  (  121,22  )  ,NP1 

121  CALL  PROBE  II  (  EL ,P I N »P0C1 ,NWB  ) 
GO  TO  2 
22    POC1  =  0.0 

2  EL  =  EL2  -  EPL 
IF(EL)28,13,13 

13  GO  TO  (  131»28  )  »NP2 

131  CALL  PROBE  II  (  EL ,P IN »POC2 ,NWB  ) 
GO  TO  8 

28  POC2  =  0.0 

8  EL  =  EL3  -  EPL 
IF(EL)25,14,14 

14  GO  TO  (  141,25  ) ,NP3 

141  CALL  PROBE  II  (  EL ,P IN »POC3 »NWB  ) 
GO  TO  5 

25  POC3  =  0.0 

5  EL  =  EL4  -  EPL 
IF(EL)26,15,15 

15  GO  TO  (  151,26  ) ,NP4 

151  CALL  PROBE  II  (  EL ,P IN ,P0C4 ,NWB  ) 
GO  TO  6 

26  P0C4  =  0.0 

6  EL  =  EL5  -  EPL 
IF(EL)27,16,16 

16  GO  TO  {  161,27  ) ,NP5 

161  CALL  PROBE  II  (  EL ,P IN ,POC5 ,NWB  ) 
GO  TO  7 

27  POC5  =  0.0 

7  EL  =  EL6  -  EPL 
IF(EL)  29,17,17 

17  GO  TO  (  171,29  )  ,NP6 

171  CALL  PROBE  II  (  EL ,P I N ,POC6 .NWB  ) 
GO  TO  30 

29  P0C6  =  0.0 
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30  CONTINUE 
GO  TO  779 

11  PNA  =  0.0 
GO  TO  882 
779  CONTINUE 

IF(  J  -  NPT 

201  POCP  =  POCN 

202  POCN  =  POCP 
800  PNA  =  POCN/ 

1  POC5  +  POC6 
882  CONTINUE 
EL  =  ELI 
CALL  PROBE 
EL  = 
CALL 
EL  = 
CALL 
EL  = 
CALL 
EL  = 
CALL 
EL  = 
CALL 


EL2 

PROBE 

EL3 

PROBE 

EL4 

PROBE 

EL5 

PROBE 

EL6 

PROBE 


)  800,201.202 


(  POC1  +  POC2  +  POC3  +  POC4  + 
) 


II  (  EL»PNA»PNL1»J 

II  (  EL,PNA»PNL2»J 

II  (  EL»PNA»PNL3»J 

II  (  EL»PNA.PNL4»J 

II  (  EL.PNA.PNL5.J 


II  (  EL»PNA,PNL6.J 


10 
881 


PN20(J)  =  PNLN/DENOMU) 
GO  TO  881 
PN20U)  =  0.0 
CONTINUE 


POC1 
POC2 
POC3 
P0C4 
POC5 
POC6 
PNL1 
PNL2 
PNL3 
PNL4 
PNL5 
PNL6 


0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 


889  CONTINUE 
RETURN 
END 

SUBROUTINE  PROD  (  EL »P I N »PNO»NWB  ) 
DIMENSION  E(102)»P(101)  »PNO(101) ,DENOM( 101) 
•   COMMON  E.  POCl,POC2»POC3»POC4»POC5»POC6» 

1  PNLl»PNL2»PNL3tPNL4,PNL5,PNL6»RN»C»A» 

2  DENOM 

DO  2  K  =1,NWB 

Z  =  E(K)  -  EL 

IF(Z)10»3»3 
3  H  =  0.3 

W  =  10.0*(Z/3.) 
6  NW  =  W  +  1.0 

Y  =  0.0 
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DO  206  L  =1»NW 

P(L)  =  C*(Z  -  Y)*PIN*EXPF<A*(Y**RN  -  Z**RN  )) 

206  Y  =  Y  +  H 
SUM  I  =  0.0 
NY  =  W 

DO  207  L. =2,  NY 

207  SUM  I  =  SUM  I  +  P(L) 

PNO(K)  =  H/2.*(P(1)  +  2.*SUM  I  +  P ( NW )  ) 
GO  TO  2 
10  PNO(K)  ■  0.0 

2  CONTINUE 
RETURN 
END 

SUBROUTINE  PROBE  II  (  EL  ,P I N ,PNO,NWB  ) 
DIMENSION  E(102),P(101) ,DENOM( 101) 
COMMON  E»POCl»POC2»POC3»POC4»POC5,POC6» 

1  PNLl»PNL2fPNL3tPNL4tPNL5tPNL6tRNtC»At 

2  DENOM 

Z  =  E(NWB  )  -  EL 
IF(Z)10,3,3 

3  H  =  0.3 

W  =  10.0*<Z/3.) 
6  NW  =  W  +  1.0 
Y  =  0.0 

DO  206  L  =  1»NW 
P(L)  =  C*(Z  -  Y)*PIN*EXPF(A*(Y**RN  -  Z**RN  )) 

206  Y  =  Y  +  H 
SUM  I  =  0.0 
NY  =  W 

DO  207  L  =  2»NY 

207  SUM  I  =  SUM  I '  +  P(L) 

PNO  =  H/2.*(P(1)  +  2«*SUM  I  +  P(NW>) 
GO  TO  2 
10  PNO  =  0.0 
2  CONTINUE 
RETURN 
END 
END 
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APPENDIX    IIC 


SUBROUTINE   GROPE 
TRIPLE   INTEGRALS 


SUBROUTINE  GROPE  (  ELI »EL2 ,EL3 » EL4, EL5 ♦ EL6 , 

1  PN1 ,PN3AN,P0CN,PNLN,EPU»EPL»EPU1,EPL1» 

2  NL1  ,NL2»NL3»NL4,NL5,NL6»NPQ,NPR, 

3  NP1,NP2,NP3»NP4,NP5»NP6  ) 
DIMENSION  E( 102) »PN1(101) ,DENOM( 101) 
COMMON  E»POCl»POC2»POC3»POC4,POC5»POC6, 

1  PNL1,PNL2»PNL3»PNL4,PNL5»PNL6»RN,C,A* 

2  DENOM 

DO  999  J  =  It  101 
IF(EU)  -  EPL)10,3,3 

3  CONTINUE 

NWC  =  10.*((E(J)  -  EPD/3.)  +  1.0 
EP  =  E(J)  -  EPL 
IF(E( J)  -  EPU) 11*4,4 

4  CONTINUE 

NWB  =  10.0*(E(J)  -  EPLD/3.  +  1.0 
IF(E( J)  -  EPU1 ) 12,5,5 

5  CONTINUE 
PIN  =  1.0 

21  EL  =  ELI  -   EPL1 
IF(EL)72,14,14 

14  GO  TO  (  141,72  ) ,NL1 

141  CALL  PROBE  II  (  EL ,P IN ,PN3A1 ,NWB  ) 
GO  TO  22 

72  PN3A1  =  0.0 

22  EL  =  EL2  -  EPL1 
IF(EL)73,15,15 

15  GO  TO  (  151,73  ),NL2 

151  CALL  PROBE  II  (  EL ,P IN ,PN3A2 ,NWB  ) 
GO  TO  23 

73  PN3A2  =  0.0 

23  EL  =  EL3  -  EPL1 
IF(EL)74,16,16 

16  GO  TO  (  161,74  ) ,NL3 

161  CALL  PROBE  II  (  EL ,P I N  ,PN3A3 »NWB  ) 
GO  TO  24 

74  PN3A3  =  0.0 

24  EL  =  EL4  -  EPL1 
IF(FL)75,17,17 

17  GO  TO  ( 171,75  ) ,NL4 

171  CALL  PROBE  II  (  EL ,P I N ,PN3A4,NWB  ) 
GO  TO  25 

75  PN3A4  =  0.0 

25  EL  =  EL5  -  EPLl 
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IF (EL) 76* 18, 18 

18  GO  TO  (  181,76  ) ,NL5 

181  CALL  PROBE  II  (  EL  ,P  I N  ,PN3A5 ,NWB  ) 
GO  TO  26 

76  PN3A5  =  0.0 

26  EL  =  EL6  -  EPL1 
IF(EL)  77,19,19 

19  GO  TO  (  191,77  )  ,NL6 

191  CALL  PROBE  II  (  EL ,P IN  »PN3A6»NWB  ) 
GO  TO  78 

77  PN3A6  =  0.0 

78  CONTINUE 
GO  TO  993 

12  PN3A  =  0.0 
GO  TO  995 
IF(  J  -  NPQ  )  993,201,202 

201  PN3AP  =  PN3AN 

202  PN3AN  =  PN3AP 
993  CONTINUE 

PN3A  =  PN3AN/(  PN3A1  +  PN3A2  +  PN3A3  +  PN3A4  + 
1  PN3A5  +  PN3A6'  ) 
995  CONTINUE 

31  EL  =  ELI  -  EPL 
IF(EL) 172,42,42 

42  GO  TO  (  421,172  ) ,NP1 

421  CALL  PROBE  II  (  EL , PN3A ,P0C1 ,NWC  ) 
CALL  PROBE  II  ( EL, 1 . 0 , POB1 ,NWC  ) 
GO  TO  32 

172  POC1  =  0.0 
POB1  =  0.0 

32  EL  =  EL2  -  EPL 
IF(EL) 173,43,43 

43  GO  TO  (  431,173  ) ,NP2 

431  CALL  PROBE  II  (  EL ,PN3A  ,POC2 ,NWC  ) 
CALL  PROBE  II  (EL,1.0,POB2,NWC  ) 
GO  TO  33 

173  POC2  =  0.0 
POB2  =  0.0 

33  EL  =  EL3  -  EPL 
IF(EL) 174,44,44 

44  GO  TO  (  441,174  ) ,NP3 

441  CALL  PROBE  II  (  EL ,PN3A,POC3 ,NWC  ) 
CALL  PROBE  II  ( EL, 1. 0 ,POB3  ,NWC  ) 
GO  TO  34 

174  POC3  =  0.0 
POB3  =  0.0 

34  EL  =  EL4  -  EPL 
IF(EL) 175,45,45 

45  GO  TO  (  451,175  ) ,NP4 

451  CALL  PROBE  II  (  EL , PN3A ,P0C4,NWC  ) 
CALL  PROBE  II  ( EL, 1. 0 ,P0B4,NWC  ) 
GO  TO  35 

175  P0C4  =  0.0 
P0B4  =  0.0 
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35 

46 
461 


176 
36 

47 
471 


EL  =  EL5  -  EPL 
IF(EL) 176,46,46 
GO  TO  (  461,176 
CALL  PROBE  I  I  ( 
CALL  PROBE  I  I 
GO  TO  36 
POC5  =  0.0 
POB5  =  0.0 
EL, =  EL6  -  EPL 
IF(EL) 177,47,47 
GO  TO  (  471,177 


)  ,NP5 

EL,PN3A,POC5,NWC 
(EL,1.0,POB5,NWC  ) 


CALL  PROBE  II 
CALL  PROBE  II 
GO  TO  48 
POC6  =  0.0 
POB6  =  0.0 
CONTINUE 
GO  TO  992 
PN1A1  =  0.0 
GO  TO  994 
IF  (  J  -  NPR 
POCP  =  POCN 
POCN  =  POCP 
CONTINUE 
PN1A1  =  POCN/M 
1  P0B4  +  POB5  + 
994  CONTINUE 
EL  =  ELI 
CALL  PROBE 


)  ,NP6 
EL,PN3A»POC6,NWC  ) 


(EL,1.0,POB6,NWC  ) 


177 


48 


11 


901 
902 
992 


)  992,901,902 


POB1 
POB6 


POB2  +  POB3  + 


EL  = 
CALL 
EL  = 
CALL 
EL  = 
CALL 
EL  = 
CALL 
EL  = 
CALL 


EL2 

PROBE 

EL3 

PROBE 

EL4 

PROBE 

EL5 

PROBE 

EL6 

PROBE 


I 


I 


I 


(  EL,PN1A1,PNL1,J  ) 


(  EL,PN1A1,PNL2,J  ) 


(  EL,PN1A1,PNL3,J  ) 


(  EL,PN1A1,PNL4,J  ) 


(  EL,PN1A1,PNL5,J  ) 


(  EL,PN1A1,PNL6,J  ) 


10 
991 


PNHJ)  =  PNLN/DENONH  J) 
GO  TO  991 
PN1 (J)  =  0.0 
CONTINUE 


POC1 
POC2 
POC3 
P0C4 
POC5 
POC6 
PNL1 
PNL2 
PNL3 
PNL4 
PNL5 


0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 


11+5 


PNL6  = 

0.0 

PN3A1 

=  0.0 

PN3A2 

=  0.0 

PN3A3 

=  0.0 

PN3A4 

=  0.0 

PN3A5 

=  0.0 

PN3A6 

=  0.0 

P0B1  = 

0.0 

POB2  = 

0.0 

POB3  = 

0.0 

POB4  = 

0.0 

POB5  = 

0.0 

POB6  = 

0.0 

999  CONTINUE 

RETURN 

END 
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APPENDIX   III 
DENSITY    FUNCTIONS 

1.  P(E)    =    C    (2   -    E)    EXP(    A(    En-    Sn)  )     ;    S 

2.  P(E)    =    C    EXP(    A    En) 

3.  P(E)    =    C 

4.  P(E)    =    C    Z    EXP(    A    (    En  -    Ethn    )) 

5.  P(E)    =    C    2    EXP(    A    En) 

6.  P(E)    =    C    2    EXP(    A    (     En  -    Sn    )  ) 

7.  P(E)    =    C    E    EXP(    A    En    ) 

8.  P(E)    =    C    E    EXP    (    A    (    En  -   Sn    )) 

9.  P(E)    =    C    E    EXP    (    A    (    En  -    Ethn    )) 

10.  P(E)    =    C/EXP(    A    En) 

11.  P(E)    =    C    E    EXP(    2(AE)n    ) 

12.  P(E)    =    C    E    EXP(    2  AE    ) 
13.P(E)    =    C    E/EXP(    AEn) 
14.  P(E)    =   C(E)3    EXP(    AEn    ) 


=    E   -Eth 
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APPENDIX    IV 
TEST 
In   order   to   test   the  program,    it   was   postulated  that 
if  all   threshold  values   were   set   equal   to   the   same  value, 
the   output   of    the   program    would  be   constant   and   the 
results    would  be   the  input    multiplied  by  that   constant. 
In   carrying   out   this   test,    the  density  function: 

P(E)  —  C    E    EXP(    A   En    ) 
was   used.      All   threshold  values   were   set    equal   to   7.2 
MEV,    the  lowest   threshold  value. 

The  result   is   a   graph   which  is   the  plot   of   the  input 
multiplied  by  a   constant   as    expected.      See   Figure  65. 
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